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ABSTRACT 

The spectral properties of type la supernovae in the ultraviolet (UV) are investigated using the 
early-time spectra of SN 200 lep and SN 2001 eh obtained using the Hubble Space Telescope 
(HST). A series of spectral models is computed with a Monte Carlo spectral synthesis code, 
and the dependence of the UV flux on the elemental abundances and the density gradient in 
the outer layers of the ejecta is tested. A large fraction of the UV flux is formed by reverse 
fluorescence scattering of photons from red to blue wavelengths. This process, combined with 
ionization shifts due to enhanced line blocking, can lead to a stronger UV flux as the iron- 
group abundance in the outer layers is increased, contrary to previous claims. 

Key words: radiative transfer — supernovae: general — supernovae: individual (SN2001ep) 
— supernovae: individual (SN 2001eh) — cosmology: observations 



1 INTRODUCTION 

Type la supernovae (SNe la) are among the most important tools to 
study the expansion history of t he Universe and t o set constraints 
on the properties of dark energy (Riess et al. 1998; Perlmutter et al. 
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Wood-Vasev et al.l 20071 : see Leib undgut 



2001; 



jPanag ia 2003: Foley et al.l 2008). Prominent examples include the 
Hubble Space Telescope (HST) spectra of SN 1992 A (5 d and 45 d 
after 6-band maximum, lKirshneretai] ri993) and the very early 
International Ultraviolet Explor er (WE) spectrum of SN 1990N 
(14d before 6-band maximum, Leibundgut et al. 1991). Recently, 
the ultraviolet/optica l telescope (UVOT) onboard the Swift satellite 



2005 for re views). With the empirical calibra- 



tion of their pe ak luminosity dPhillip j[l993l ; IPhillips et all 1 19991 ; 
|jha et alj|2007l) , it is feasible to probe the properties of the Uni- 
verse out to a redshift of; » 1 and highe r ( Gold haberetai]|200ll ; 
iTonrv et alfeoOSURiess et al .11 20041. 1 20071) . Understanding the for- 
mation of the ultraviolet (UV) spectrum of SNe la is important for 
two reasons: the optical spectrum is strongly affected by line block- 
ing in the UV, and optical observations of high-redshift supernovae 
primarily cover the rest-frame UV bands. Comparisons between lo- 
cal and distant SNe la have recently identified the UV as an impor- 
tant region to search for possible effects o f evolution dFolev et "all 
l2007l ; lBronder et al]2008f ; lEliis et alj|2008l) . 

Efforts have been made to increase the sample of local SNe la 
observed in the UV by space-borne telescopes. The amount of data 
available for SNe la in the UV is very small and mostly limited 
to single epochs - usually around (or well after) maximum light 
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Gehrels et al. 2004) has also provided valuable UV data on SNe la 



Immler et al. 2006). 



To interpret such data correctly, it is necessary to understand 
the formation of the UV flux in SNe la and its sensitivity to physical 
parameters of the explosion and ambient medium. Several groups 
have studied the properties of t he UV flux using radiative transfer 



have studied the properties ot the U V tlux using radiative transfer 
models of differing complexity jPauldrach et al.ll 1 996 ; Fish er et all 
ll997l;lHoflich et aljl9 98: Mazzali 200Ct lLentz et alfeOOll) . 



In this paper, we use radiative transfer models to study the 
properties of SN 2001 eh and SN 200 lep, which were observed as 
part of the Supernova INtensive Study (SINS; GO-9114; PI R. P. 
Kirshner) with HST. Based on the spectral models, we investigate 
how a change of physical parameters may affect the observed UV 
flux of a SN la. In particular, when using SNe la as standard can- 
dles to measure cosmological parameters, it is important to know 
if parameters that are likely to change over a large range of red- 
shifts can introduce systematic errors in the distance measurement. 
Among the parameters that may be affected by the environment of 
SN la progenitors is the abundance of heavy elements in the outer 
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shells of the ejecta. These layers are of particular importance for 
the formation of the UV spectrum during the photospheric phase. 
The two supernovae chosen for this study are at the opposite ends 
of the luminosity distribution of normal SNe la. 

From spectral analysis of SNe la, we know that the compo- 
sition is dominated by a mixture of intermediate-mass elements, 
such as Si, S, and Mg, and iron-group elements (Ni, Co, and Fe). 
At early epochs, until a few weeks after maximum brightness, the 
spectrum of a SN la consists of a number of broad absorption 
troughs that generally originate from blends of v arious lines. (For 
a re view of the observ ed properties of SNe la, see Filip penko|[l997l 
and Leibundgut 2000.) The emission of SNe la is entirely powered 
by the radioactive decay of 56 Ni created in the explosion. The y-ray 
photons from radioactive decay undergo Compton scattering, cre- 
ating fast electro ns which deposit their energy via collisions in the 
expanding ejecta dTruran et al.ll967l ; IColgate & McKeel 19691) . The 
non-thermal continuum seen at early epochs around maximum light 
is formed by the overlap of thousands of weaker lines. The spec- 
trum peaks at a wavelength between 4000 and 5000 A and drops 
steeply in the UV, bluewards of the Ca H&K lines, which typically 
have an absorption m inimum at about 3750 A near maximum light 
jBlondin et alj2006h . 

The UV wavelength bands of SNe la comprise only a small 
fraction of the total emitted light, and are generally characterized by 
the absence of strong features. Singly and doubly ionized species of 
heavy elements such as Ti, Cr, Fe, Co, and Ni have a large number 
of line transitions in the UV bands that effec t ively block most of 
the emission in these bands dKarp et alJI 19771; iHoflich et alJI 19931 ; 
IPauldrach et alJI 19961) . IlucvI dl999bh and lMazzalil i feOOOl) find that 
essentially all UV photons created at deeper layers in the ejecta 
are immediately absorbed and down-scattered to lower energies 
(longer wavelengths) where the line opacity is smaller. The cru- 
cial mechanism for the emission of photons in the UV is found to 
be a reverse-fluorescence process where absorption o f red photons 
is followed by emission in the blue and UV bands dLucvf l 1999b; 
Mazzali 2000). These processes have a significantly lower proba- 
bility of occurrence but in the outer region of the ejecta the emission 
of bluer photons via reverse fluorescence can be significant because 
there are many more red photons available than blue ones such that 
the lines in the blue are not saturated. Blue photons can escape 
from a line rather than getting re-absorbed locally. The UV spec- 
trum is predominantly shaped by this fluorescence spectrum orig- 
inating from a confined layer of high-velocity material, which is 
small compared to the dimensions of the expanding ejecta dMazzalil 
2000). The shape of the emission spectrum is determined by the 
atomic properties of the ions in the emission layer. This makes the 
UV an interesting wavelength band in which to study possible pro- 
genitor scenarios because the progenitor star may leave a charac- 
teristic imprint on the composition of the outermost layers of the 
exp anding shell. 

iLentz et all d200ll) investigated the dependence of the UV flux 
of SNe la on the composition of the outermost layers of the ejecta, 
adopting different fractions of solar metallicity for the material 
above v ~ 14 000 km s~'. They found that overall the UV flux de- 
creases with increasing metallicity. They also pointed out, however, 
that the larger number of strong absorption lines with higher metal- 
licity can cause a back reaction of the blocked radiation field on 
the temperature and ionization structure (line blanketing), which is 
difficult to disentangle from the line-blocking effect. Therefore, the 
UV flux may vary non-linearly with a change of composition. 

In this paper we investigate the sensitivity of the UV spec- 
trum to modifications of the composition in the outer layers. We 



also consider variations in the density structure of the underlying 
explosion model. 

The paper is organized as follows. In Section [2] we present 
observations of SN 200 lep and SN 2001 eh covering the UV wave- 
length range as an example of SNe la with different UV properties. 
Section[3]outlines the radiative transfer code used for this work and 
an analysis of the observed spectra with the help of synthetic spec- 
tra. The radiative transfer models for both supernovae are used as a 
basis for a series of models to study the properties of the UV flux 
with respect to changes in composition and density structure. The 
results of this series and the physical processes responsible for the 
formation of the UV flux are discussed in Section|4] Section[5]gives 
the conclusions of our study. 



2 OBSERVATIONS 

We present observations and synthetic spectra of SN 2001 eh and 
SN 200 lep as examples of SN la that exhibit significantly different 
spectral properties, especially in the UV. A more detailed discus- 
sion of the complete data available for the two supernovae (UBVRI 
light curves, optical spectra, and one additional UV spectrum for 
each SN) will be presented in a forthcoming paper (Blondin et al., 
in preparation). Here we present the data relevant to the study of 
the UV spectra discussed in this paper. The observations are sum- 
marized in Table Q] 

SN 2001eh was discovered on 9 Sep. 2001 (UT dates are 
used throughout this paper) by M. Armstrong dHurst & Armstrong! 
1200 lh and independently b y M. Ganeshalingam and W. D. Li 
dGaneshalingam et al.ll200l|) with th e Katzman Automatic Imaging 
Telescope (KAIT; iFilippenko et alj|200lh as part of the Lick Ob- 
servatory and Tenagra Observatory Supernova Search (LOTOSS). 
It occurred in UGC 1162, which have a NEEQ recession veloc- 
ity of cz = 1113 1 kms~' (from The Updated Zwicky Catalog; 
iFalco et al.l 1 1999h . A CCD spectrum obtained on 11 Sep. with 
the Shane 3-m reflector at Lick Observatory revealed that SN 
2001 eh was a SN la about a week before maximum brightness 
(Ganeshalin gam et alfeOOll) . 

SN 2001ep was discovered on 3 Oct. 2001 bv lHutchings & Lil 

d200ll) as part of LOTOSS. It occurred in N GC 1699, which has a 
NED recession velocity of cz = 3901 km s -1 (de Vauc ouleurs et al.l 
Il99lh . A CCD spectrum obtained on 10 Oct. with the Fred L. 
Whipple Observator y (FLWO) 1.5-m Til linghast telescope and 
FAST spectrograph dFabricant et al.lll998|) showed th at it was a 
SN la near maximum brightness ( Matheson et al. 200j]). 

SN 2001eh and SN 2001ep were observed with HST (Cycle 
10) as targets of opportunity. The data were obtained with the Space 
Telescope Imaging Spectrograph (STIS) using a 52"x0.2"slit. For 
both supernovae we used the G430L grating to obtain "optical" 
spectra (2900-5700 A), the G230L grating to obtain near-UV spec- 
tra (1600-3100 A), and the G140L grating to obtain far-UV spectra 
(1150-1700 A). No clear signal was detected for any observation 
taken with the latter setup. The reductions followed the standard 
STIS reduction pipeline. Spectra were extracted from the processed 
images using the IRAf| SPECRED long-slit reduction packet. 



' The NASA/IPAC Extragalactic Database (NED) is available at 
http://nedwww.ipac.caltech.edu 

2 http: //tdc-www. harvard. edu/uzc/ 

3 IRAF is distributed by the National Optical Astronomy Observatories, 
operated by the Association of Universities for Research in Astronomy, Inc., 
under contract to the National Science Foundation of the United States. 
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Figure 1. KAIT B-band light curves of SN 2001eh and SN 2001ep. The in- 
set shows both light curves out to ~20 d past B-band maximum, normalized 
in magnitude at the peak. The dotted line at +15 d highlights the difference 
in Ami5(B) between the two supernovae (see text for details). 



The optical spectra presented in this paper were obtained with 
the FLWO 1.5-m FAST spectrograph using a slit of width 3". For 
both supernovae we used a 300 line mur 1 grating with two tilts: 
one covering -3700 - 7500 A, the other -5200 - -9000 A. The 
spectra were reduced and calibrated employing standard techniques 
in IRAF and our own IDL routines for flux calibration. 

The KAIT B-band light curves for SN 2001 eh and SN 200 lep 
are shown in Fig. [T] Both are sufficiently well sampled to enable 
a direct measurement of the decline-rate parameter Amis(B). For 
SN 2001eh we find Am 15 (B) = 0.71 ± 0.05, while for SN 2001ep 
we find Am 15 (B) = 1.41 ± 0.04. Thus, SN 2001eh is a "slow de- 
cliner," while SN 2001ep is at the fast end of typical "intermediate- 
decliner" SNe la. 

The differences in light-curve shape between the two super- 
novae are reflected in their spectra. Using the Supernova I dentifi - 
cation (SNID) cross-correlation code of iBlondin & Tonrvl j2007t) . 
we find that the best-match for SN 2001 eh is the overlum inous 
SN 1999aa (Am ls (B) = 0.75 ± 0.02; iKrisciunas et alj|2000l) , and 
that the best-match template spectrum for S N 2001ep is the "tran - 
sitional"SN2004eo(Am 15 (B) = 1.45±0.04: |Pastorello et alj2007h . 

We fitted the complete UBVRI photom etric data with the 
MLCS2k2 light-curve fitter o flJha et alj J2007h . The time of B-band 
maximum is found to be 2001 Sep. 17.53 (JD = 2452170.03±0.71) 
for SN 2001eh and Oct. 17.72 (JD = 2452200.22 ± 0.45) for 
SN 200 lep. The derived line-of-sight extinction in the host galaxy 
is A v = 0.08 ± 0.06 and A v = 0.43 + 0.12 mag, respectively. 
The Milky Way foreground reddening is small for both supernovae 
{E(B - V) = 0.06 mag for SN 20 leh and E(B -V) = 0.05 mag for 
SN 2001ep; ISchlegel et al]l998T) . 

The UV and optical spectra of each SN were combined into 
a single spectrum. The far-UV data taken with the G140L grating 
were not included in either of the combined spectra given the lack 
of signal. The mean UT date for the SN 2001eh spectrum is 2001- 
09-25.93, corresponding to -9 d past B-band maximum. The mean 
UT date for the SN 2001ep spectrum is 2001-10-26.88, also corre- 
sponding to -9 d past B-band max i mum. Assuming a B-band rise 
time of -19.5 d dConlev et al.l2006l ; lRiess et aljl999l) . both spectra 



correspond to 28-29 d past explosion. In what follows we assume 
that both spectra are at t = 29 d past explosion. 



3 SPECTRAL MODELS 

3.1 The radiative transfer model 

To model synthetic spectra of SNe la at early epochs, we use a 
Monte Carlo spectral synthesis code with an approximate descrip- 
tion of non-LTE (i.e., not assuming local thermodynamic equi- 
librium, LTE). Here we outline the most important aspects of 
our method as they are relevant for the d iscussion in thi s work . 
More detailed descripti ons c an be found in Abbot t & Lucvl dl985l), 
IMazzali & Lucvl dl993h . and lLucvl dl999bl) . as well as in lMazzalil 
d2000h . This last reference also gives the sources of atomic data 
used by the code. 

For the underlying explosion model we take the den- 
sity structure from the one-dimensional deflagration model W7 
dNomotoetalJl984h . The composition is adjusted to match the ob- 
served spectrum. For this study we use the version of the code that 
allows a r adial variation of the composition within the ejecta de- 
scribed by Steh le et al.1 d2005h . 

The code is based on the Schuster-Schwarzschild approxima- 
tion - that is, we assume that the radiative energy is entirely emitted 
at the inner boundary of the computational grid with the spectral 
shape of a blackbody. The spectra are calculated assuming station- 
arity and ignoring any time-dependences. This is a good approxi- 
mation at l east for ID mode ls of SN la later than a few days after 
explosion dKasen et alj|20061) . We do not consider the deposition 
of energy by y-ray photons from the 56 Ni and 56 Co decays within 
the atmosphere, nor do we try to derive the total luminosity di- 
rectly from the 56 Ni mass. Therefore, luminosity is a free param- 
eter. To constrain the model, the location of the inner boundary 
("photosphere'0) has to be specified in terms of a velocity vn, and 
the elapsed time t since the explosion, assuming homologous ex- 
pansion (r(t) = vt) of the ejecta. In addition to those physical input 
parameters, observational parameters such as distance and redden- 
ing have to be specified to allow direct comparison of the synthetic 
and the observed spectrum. 

The code uses an approximate description of non-LTE to de- 
rive the atomic level populations and the ionization equilibrium. 
Collisional and continuum processes are not considered, apart from 
electron scattering. Since the opacity in SN la ejecta is largely dom- 
inated by the opacity of Doppler-shifted and broadened spectral 
lines, the true continuum plays only a minor role. Given the low 
densities in the ejecta, collisional processes are also not expected 
to have a major impact on the state of the gas. It is not possible 
to derive a consistent gas tempera ture in this approximation. Thus, 
following IMazzali & Lucvl d 1993b . the local gas temperature is ap- 
proximated to be 90% of the radiation temperature Tr of the radia- 
tion field, which is determined from the mean energy x of a photon 
in a Planckian radiation field, 



4 We generally avoid the term "photosphere" for this inner boundary. In 
contrast to stars, where this term refers to the radius at which the atmo- 
sphere becomes optically thick with respect to a wavelength-independent 
mean optical depth, this concept cannot be used in hydrogen-deficient su- 
pernovae, where the optical-depth scale is strongly wavelength-dependent 
owing to the shallo w density gr adient and the absence of strong continuum 
opacities (see, e.g. JSauer et alj2 006). 
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Table 1. UV and optical spectra of SN 2001eh and SN 2001ep 



IAU name 


UT Date 


Telescope+ 


Slit 


Grating 


Range 


Exptime 






Instrument 


(") 




(A) 


(s) 


2001 eh 


2001-09-25.67 


HST +STIS 


52x0.2 


G140L 


1140-1730 


1340 


2001 eh 


2001-09-25.76 


HS 7/+STIS 


52x0.2 


G230L 


1570-3180 


2x2700 


2001 eh 


2001-09-25.66 


HS 7/+STIS 


52x0.2 


G430L 


2900-5700 


500 


2001 eh 


2001-09-26.35 


FLWO 1.5m+FAST 


3 


300 


3720-7540 


1200 


2001eh 


2001-09-26.37 


FLWO 1.5m+FAST 


3 


300 


5200-9040 


1200 


200 lep 


2001-10-27.99 


HST +STIS 


52x0.2 


G140L 


1140-1730 


1630,2690 


2001ep 


2001-10-28.09 


HS 7+STIS 


52x0.2 


G230L 


1570-3180 


2690 


200 lep 


2001-10-27.95 


HS T+STIS 


52x0.2 


G430L 


2900-5700 


200 


200 lep 


2001-10-25.42 


FLWO 1.5m+FAST 


3 


300 


3720-7540 


1200 


2001ep 


2001-10-25.44 


FLWO 1.5m+FAST 


3 


300 


5200-9040 


1200 



liv 



£° vB r dv 360 



k B T kuT fB v dv x 4 



£(5) a 3.832 



(1) 



(see also lLucvlll999ah . The temperature of the blackbody at the 
inner boundary is determined in an iteration cycle, using the con- 
straint that the desired luminosity L is obtained at the outer bound- 
ary. Iteration is necessary because the amount of radiation that is 
scattered back into the core is not known a priori. 

A crucial quantity for the derivation of the ionization and exci- 
tat ion equilibrium is the "e quivalent dilution factor" W introduced 
bv lMazzali & Lucvl ( fl993h . This factor, which is derived from the 
nebular approximation for the radiation field in an optically thin 
medium, is used to modify the LTE Saha-Boltzmann excitation for- 
mula. While in a gaseous nebula W describes a purely geometrical 
dilution of the radiation field, the procedure developed for super- 
nova ejecta derives W directly from the assumed relationship be- 
tween the Planck function at the radiation temperature T R and the 
radiation field /: 



J = WB(T) = W-TL 

71 



(2) 



where a denotes the Stefan-Boltzmann constant. This approxi- 
mation has been shown to produce a reasonably good represen- 
tation of non-LTE ion ization and excitation in supernova ejecta 
dMazzali&Lucvll 19931) . 

Line process es are treated usi ng the lin e-branching formal- 
ism introduced by Lucv ( 1999b) and Mazzali (2000). While earlier 
work with this code assumed pure resonant scattering of photons 
in spectral-line transitions, this improvement allows photons to be 
re-emitte d in a diff e rent tr ansition than the one by which they were 
absorbed. iMazzalil q2000f) showed that the inclusion of branching 
processes is necessary to reprod uce the flux in the UV. 

Following the description o f lLucvUl999bl) , the emergent spec- 
trum is derived by a formal solution of the transfer equation using 
line and continuum source functions that are extracted from the 
Monte Carlo simulation. The advantage of this method over using 
emerging Monte Carlo packets directly is a significantly reduced 
number of Monte Carlo packets to obtain a converged spectrum 
with low Monte Carlo noise. Because the UV flux of a SN la is in- 
trinsically very low, one must nevertheless ensure that the number 
of packets used in the Monte Carlo experiment is sufficient to pro- 
vide a converged solution for the source functions. We performed 
convergence tests in the UV and determined that the relative scat- 
ter of the spectrum caused by Monte Carlo noise in the wavelength 
range between 1500 and 3000 A is generally well below 10% for 



500 000 photon packets. With this number the scatter in the optical 
spectrum is less than 0.4%. 



3.2 Synthetic spectra for SN 2001ep and SN 2001eh 

The models and the observed spectrum of SN 2001 eh are shown in 
Fig.|2ji; those of SN 2001ep are shown in Fig. |2j). The inserts in 
the figures show the UV part of the spectra on a logarithmic scale. 
The direct comparison of both observed spectra clearly shows that 
these supernovae exhibit fairly different spectral appearances. In 
SN 2001eh the UV flux blue of Ca II H&K at -3750 A is stronger 
with relation to the optical spectrum. The characteristic line fea- 
tures are fairly narrow and are found at lower velocities than in the 
spectrum of SN 200 lep. 

Modelling the spectra requires knowledge of the distance and 
reddening of the supernovae. For the distance we use the data avail- 
able for the respective host galaxies listed in NED. The redshift 
of the host galaxy of SN 2001eh (UGC 1162, z = 0.037) cor- 
responds to a distance mo dulus of m - M = 35.9 mag assuming 
H = 73kms 1 Mpc 1 dRiess et alj|2005t) . NGC 1699, the host 
galaxy of SN 2001ep, has a redshift of z = 0.013, resulting in a dis- 
tance modulus of m - M = 33.6 mag. For the model of SN 2001eh 
we assume a total reddening of E(B - V) = 0.07 mag, while for 
SN 2001ep we use E(B-V) = 0.15 mag. The value for the red- 
dening of SN 2001ep is at the lower end of the uncertainty range 
estimated from observational indicators (see Section^; higher val- 
ues for E(B-V), however, require a higher luminosity and lead to a 
shift in the ionization towards doubly ionized species of Fe and Co 
that are not seen in the observed spectrum. 

For both supernovae the epoch was taken to be 29 d after 
explosion (see Section O. The synthetic spectrum of SN 2001eh 
was modelled with a luminosity of L = 1.15 x 10 43 ergs~'. For 
the spectrum of SN 2001ep we need a lower luminosity of L = 
7.98xl0 42 ergs -1 . For both mode ls we adopted the density structure 
from the W7 explosion model dNomoto et al.ll 1984h : the composi- 
tion was adjusted to fit the observed spectrum. To model SN 2001eh 
we set the velocity of the inner boundary to v ib = 5200 km s . For 
SN 2001ep we used v ib = 5000 km s~'. Above the inner boundary, 
the ejecta were divided into several radial zones where the com- 
position was adjusted independently. The composition of the outer 
zones was determined using the spectra of earlier epochs available 
for both supernovae. The outcome of theoretical explosion models 
such as W7 and the results of tomography studi es of SN la with 
similar propertie s (SN 2002bo Istehle et alj|200lh for SN 2001 eh 
and SN 2004eo dMazzali et alj|2008h for SN 200 lep) served as a 
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Figure 2. Synthetic and observed spectra of (a) SN 2001eh and (b) 
SN 200 lep, both at t = 29 d after explosion. The inserts show the UV part of 
the respective spectra with a logarithmic scale. The labels indicate the ions 
that dominate the formation of the various absorption features. The blue and 
UV parts of the spectrum are formed by a blend of a large number of lines, 
precluding a complete identification of the contributors. 



rough guideline. The inner region (v < 8000 km s~ 1 ) of the ejecta of 
SN 2001eh contain mostly (> 80%) elements generated by burning 
to nuclear statistical equilibrium (NSE; Fe-group material, mostly 
56 Ni). This is followed by a layer containing ~60% intermediate 
mass elements (IME) (Si, S, Mg, Ca) and -20% of NSE material 
out to 1 2 000 km s~ 1 . The layers above 1 2 000 km s~ 1 are dominated 
by unburned O and C (58-98%) with small amounts of IME (be- 
tween 2 and 20%) and NSE material (0.5 -4.5%), decreasing out- 
wards. In the less luminous SN 2001ep the IME-rich region extends 
to lower velocities (-6000 km s~') and the O-rich layer reaches 



down to -9000 km s" 1 . The mass fraction of NSE material appears 
to drop rapidly (< 25%) above ~7000 km s _1 . In the outermost lay- 
ers above -10 000 km s" 1 the mass fraction of IME material drops 
from -40 to 1% outwards, NSE material from 3 to 0.2%; The NSE 
abundance in the outer layers is so mewhat lower than in SN 2001eh 
which is in line with the findings o f iMazzali et at] (2007 ) . The pres- 
ence of small amounts of NSE material in the outermost layers is, 
however, necessary to fit the UV part of the spectrum in both su- 
pernovae. The original W7 density structure only extends to a ve- 
locity of 24 000 km s~'. To ensure that the choice of the maximum 
velocity of the computational grid does not affect the outcome of 
the model, the density structure is extrapolated to higher velocities 
using the power-law index of the outer zones of the input-density 
structure. In the model series discussed later in this paper, the maxi- 
mum velocities is set to large values (-70 000 km s _I ) to avoid arte- 
facts generated by this boundary. 

The models fit the spectra fairly well in the UV and optical 
parts of the spectrum. The offset in the wavebands redwards of the 
Si II feature at -6100A is caused by the breakdown of our as- 
sumption of a thermal lower boundary. The opacity at these wave- 
lengths is dominated by electron scattering, which does not mod- 
ify the spectral properties of the radiation field. The radiation at 
these wavelengths originates from layers below our adopted inner 
bound ary and is sh aped by a large number of individually weak 
lines dKasenl 120060 . With our assumption of stationarity and the 
chosen set of atomic data, the shape of the pseudo-continuum can- 
not be correctly describe d by the model even if extended to deeper 
layers dSauer et"al]|2006l> . Additionally, a description of the ejecta 
below the lower boundary implies that we rely on specific explo- 
sion models, which introduce additional uncertainties. The absorp- 
tion features blue of -2700 A in the spectrum of SN 200 lep (in- 
sert of Fig. [2p) appear to be too blue in the synthetic fit. This mis- 
match probably indicates that the outer part of the density structure 
used to compute the model is not appropriate. This is not surprising 
given that the W7 model generated more 56 Ni than what is expected 
for SN 200 lep and may therefore also have a larger kinetic energy 
than required to describe this SN. Another indication of an incorrect 
density structure is that the model fails to reproduce the separated 
Fe II features around 5000 A which produce the apparent emission 
feature at -4950 A. 

The major absorption features are labelled in Fig.[2] However, 
especially towards the blue and UV, where the features are blends 
of many different lines, ions other than those indicated in the plot 
also contribute to the shape of the absorption. SN 2001 eh is bluer 
than SN 200 lep and has a different spectral appearance at the same 
epoch. The spectrum of SN 2001 eh contains a significant contribu- 
tion from doubly ionized species such as Fe III and Co III, while 
the spectrum of SN 2001ep is dominated by lines of singly ionized 
species. The characteristic Fe II features (e.g., between 4500 A and 
5000 A) are significantly weaker in SN 2001eh than in SN 2001ep. 
The comparison of the dominant ionization stages indicates that 
the gas temperatures in the ejecta of SN 2001eh are higher than in 
SN 2001ep. In the models of SN 2001ep the temperatures at simi- 
lar velocities are generally 500-1000 K lower than in the SN 2001 eh 
model. 

More differences are visible in the region bluewards of Ca II 
H&K. Relative to the optical emission peaks, SN 2001eh has sub- 
stantially more flux in this region than SN 2001ep. Further to the 
UV the Mg II 2800 A absorption, which dominates the shape of the 
spectrum in the region around 2600 A, is stronger in SN 2001ep. 

The differences seen in the UV result from different composi- 
tions or densities of the outermost layers of the ejecta. In particular, 
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Fe-group ions from Ti to Ni have a large number of spectral lines in 
the UV that lead to large optical depths out to the highest velocities 
where this part of the spectrum is shaped. 



4 THE UV FLUX IN SN la — A MODEL SERIES 

using a series of parameterized models where we modified first 
the composition and then the density gradient in the outer layers. 
Composition is expected to be the prima ry quantity respon sible for 
changes in the UV spectrum (see e.g., lLentzetal.ll200lh . In ad- 
dition, the imprint of a progenitor is likely to manifest itself as a 
variation of the composition in the outer layers. A modified den- 
sity structure can also be responsible for changes in the appearance 
of the UV spectrum. Different explosion models predict variations 
in density, in particular at the higher velocities. Specifically, mod- 
els with a deflagration-to-detonation transition (DDT) occurring at 
some stage during the explosion accelerate the matter to higher ve- 
locities resulting in higher d ensities and more burned material in 
the outer layers of the ejecta dKhokhlovl 1 19911; iHoflich et al.ll 19951 ; 
llwamoto et ail 19991 ; iRopke & Niemeverj|2007l ; iRopkell 20071) . 

4.1 Variation of the composition 

4.1.1 Setup 

We use the models discussed in the previous section to set up a pa- 
rameterized model series in which the abundance of heavy elements 
is varied stepwise in the outermost zone. For the first test series we 
keep the density structure of the model unchanged and balance the 
modification of heavy-element abundances with a change of the 
oxygen abundance. In the outer layers of the models, O is the most 
abundant element and causes very few visible features in the early- 
time spectrum of SNe la. Therefore, even drastic changes in the 
oxygen abundance do not affect the spectral appearance; the mod- 
els primarily reflect the effects of a variation in the abundance of 
heavy elements. (For the most extreme cases a change of the elec- 
tron density may become significant. A composition dominated by 
oxygen can provide more electrons per unit mass than IME or Fe- 
group elements. Most of this effect is, however, counterbalanced 
by the higher ionization potential of oxygen resulting in a larger 
fraction of neutral atoms than for heavier elements.) 

Within each series all model parameters are kept constant and 
only the abundance of a specific element or group of elements is 
increased or decreased above a velocity u S h c ii- For the models of 
SN 2001eh we set u shcll = 15 000 km s _1 , while for the SN 2001ep- 
based models we use u s h e ii = 14 500 km s -1 . The velocities are cho- 
sen to be roughly the range at which the radiation field in the op- 
tical wavelength bands is mostly decoupled from the ejecta mat- 
ter (cf. Fig [6j». The compositions of the outer layers for the two 
base models are shown in Table [2] The models are calculated out 
to v ~ 70 000 km s~ where the density is low enough to ensure that 
the outer cutoff velocity does not affect the model results. As for the 
original models, we use 500 000 Monte Carlo packets to keep the 
Monte Carlo-induced noise in the UV spectrum sufficiently low. 

We chose to vary the composition of (i) Fe alone (correspond- 
ing to stable Fe not produced via the S6 Ni-decay chain), (ii) 56 Ni 
(at the epoch of 29 d this corresponds to 3.7%Ni 79.7% Co and 
16.6% Fqfj), and (iii) Ti and Cr together. The last two elements do 
not produce strong visible features in the spectrum, but affect the 

5 using half lives from Junde ( 1999) 



Table 2. Outer-layer composition (mass fractions) of the two base models. 
Fe listed here corresponds to Fe not generated in the 56 Ni-decay. The decay 
products of 56 Ni after t = 29 d are listed separately below. 
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UV flux because they have a large number of relatively weak lines 
that densely cover the UV and therefore contribute significantly to 
the total opacity. For each base model we increased and decreased 
the composition of the groups of elements in steps of 0.5 dex. 

4.1.2 Results 

The spectra of the series are shown in Fig. [3] The models in the left 
column are based on SN 2001eh, the ones in the right column on 
SN 2001ep. Each column shows the model for Fe variation alone 
in the upper plot, the 56 Ni series in the middle, and the Ti/Cr series 
in the bottom plot. The insets show the UV part of the spectrum in 
logarithmic scaling for each set. The colour coding represents the 
different models - the shades of blue represent models in which 
the abundance of the element groups have been decreased, while 
the shades of red are models in which the abundances have been 
increased. For comparison the original model is shown in black. 

Overall, an increase of abundances has a stronger effect than 
a decrease. This is not surprising, because the abundance of heavy 
elements was already relatively low in the outer layers of the base 
models. The comparison of the three groups of models shows that 
the relative change of the optical part of the spectrum is small to 
moderate; the variation of Fe alone has the strongest effect on the 
optical part especially in the series of spectra based on SN 200 lep 
(Fig. [3] right panels). 

The variation of the blue and the UV part of the spectra dif- 
fers among the different groups of models. In the SN 2001eh series 
the absorption feature at 3050 A becomes deeper with higher Fe 
abundance; in the SN 2001ep series this feature does not change 
significantly. The region red of this feature, up to ~3500 A, shows 
significantly enhanced flux for higher Fe abundance in both cases. 
In the SN 200 lep series this higher flux level also affects the re- 
emission peak bluewards of the Ca II H&K feature, which remains 
mostly unchanged in the SN 2001eh models. The models where Fe 
is decreased are very close to the base model in both series. 

Further to the UV, bluewards of ~ 2500 A, the models based 
on SN 2001eh (Fig. [3] upper-left panel) show little change with the 
variation of the Fe abundance. In the corresponding SN 200 lep- 
series (Fig. [3] upper-right panel), we observe the trend that models 
with higher Fe abundance have more flux in the UV, although the 
effect here is also small. 
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Figure 3. Model series for varying composition of the outer layers for the SN 2001eh-based model series (left column) and the SN 2001ep-based model series 
(right column). The upper panels show the series where Fe alone is varied, the middle panels show the series with a variation of 56 Ni (i.e. Ni, Co, and Fe 
together), and the lower panels show the series for Ti and Cr. Colour coding indicates different models — blue lines represent models where the abundance of 
the respective element is decreased, red lines show models where the abundance is increased. The black line refers to the original model. 



The variation of 56 Ni is visible in a similar region as in the Fe 
models, but with the opposite trend. Models with more 56 Ni have a 
deeper absorption at ~3300 A, which is mostly caused by stronger 
Co III lines. Consequently, the emission peak bluewards of the Ca II 
absorption becomes stronger in those models (this peak is largely 
induced by the P-Cygni re-emission of those Co-lines). Again the 
variation in this peak is stronger in the SN 2001ep-based models. In 
the UV part blue of about 2500 A the effect of decreasing the 56 Ni 
abundance leads to more flux while increasing the 56 Ni abundance 
leads to less flux. 

The bottom panels of Fig.[3]show the models where the abun- 
dances of Ti and Cr have both been varied. The most prominent re- 
sult of this variation is the change of the double-peaked re-emission 
feature centred at 3000 A. For low Ti and Cr abundances in the 
outer layers this double-peaked structure becomes more and more 
prominent, while it almost entirely disappears for high mass frac- 
tions of Ti and Cr. This behaviour is more pronounced in the spectra 



of the SN 2001 eh- series. In both series high Ti and Cr abundances 
result in an additional absorption feature at 3050 A, mostly caused 
by Ti II. Further to the UV the trend again reverses, leading to 
more flux for the models with higher Ti and Cr abundance. Blue 
of ~2700 A reductions of Ti and Cr have only a small effect. In the 
optical spectrum Ti and Cr predominantly affect the blue wing of 
the deep absorption around 4300 A. 

To quantify the change in the models better we integrate the 
flux in different bandpa ss filters of the WFPC2 instrument and the 
standard (/-band filter jBesselllll990h . Fig. H shows the logarith- 
mic representation of the Fe-model series of the SN 2001ep-based 
models in Fig. [3] (upper-right panel) together with the shape and 
position of the different bandpass filters in the upper panel. 

Fig.[5]shows the integrated flux of the model spectra as a func- 



See http://www. stsci . edu/hst/wfpc2 
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Figure 4. Logarithmic plot of the SN 200 lep model series with varying 
Fe content (cf. Fig. [3] upper-right panel). The upper panel shows the posi- 
tion and shape of the WFPC2 bandpass filters used for the comparison of 
integrated flux shown in Fig|5] 



tion of the change in the abundance of the respective set of elements 
relative to the base model. The different lines and symbols corre- 
spond to the filters indicated in the legend. Smaller values of Am 
correspond to more flux in the corresponding filter. Overall, the 
model series for the two base models show similar behaviour, al- 
though the absolute numbers are slightly different. The plots con- 
firm what is visible in the spectra: the increase of Fe and Ti/Cr 
leads to a higher flux in some passbands, while the modification 
of 56 Ni has the opposite effect. The U band is almost unaffected 
by a decrease of heavy-element abundances, while an increase al- 
ways leads to more flux (although to a small degree in some cases). 
The feature that most strongly affects the £/-band is the re-emission 
peak bluewards of the Ca II H&K absorption. The largest variation 
is usually seen in the filter centred around 1850 A (f!85w). 



4.1.3 Discussion 

The change in the UV flux with varying composition is mediated 
by two main effects. Firstly, an increased content of heavy elements 
with a large number of lines in the UV increases the probability 
for the reverse fluorescence to occur. Secondly, the increased line 
opacity leads to a stronger line blocking effect, which affects the 
temperature structure and consequently the ionization structure in 
the ejecta. 

The two plots in Fig.|6]illustrate the reverse fluorescence effect 
in the Monte Carlo radiative transfer model for SN 2001ep. The 
points indicate the initial wavelength at which a photon package has 
been created versus the final wavelength at which it emerges from 
the envelope after the last scattering event. The grey scale indicates 
the velocity of the shell where the last interaction occurred. Darker 
points correspond to higher velocities. Resonance scatterings that 
do not change the wavelength of the package are located on the 
diagonal. 

The data in the upper panel are from the simulation of the 
base model for SN 2001ep; the data in the lower panel are from 
the model with high Fe abundance in the outer shells. Both simula- 
tions have a similar total number of emerging packets (~ 180 000), 
but the diagrams show only a fraction of this total number because 
many packets are located on identical positions in the wavelength 
grid. In the model with higher Fe abundance in the outer layers the 
photons encounter more line events and are therefore more likely 
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Figure 6. Fin al versus initial wavelength of all photon packets that emerge 
to infinity (cf. Mazzali 2000; Pinto & Eastman 2000). The grey shading in- 
dicates the velocity of the last scattering event. The darker a point is, the 
higher the velocity of the shell at which the photon packet encountered 
the last interaction with a line. (For wavelength pairs where more than one 
photon packet escapes, only the one with the highest velocity is shown. 
Note that a single point usually represents a larger number of photon pack- 
ets.) The diagonal line represents all photons that are resonance scattered. 
All points above the diagonal are packets scattered to redder wavelengths, 
while all points below the diagonal represent packets that are scattered to 
bluer wavelengths through the reverse-fluorescence process. The upper plot 
shows the base model of SN 2001ep; the lower plot shows the same quanti- 
ties, however, for the model where the Fe abundance in the outer shell above 
14 500 km s was enhanced by 2 dex. The overall comparison shows that in 
the latter case, more photon interactions occur in the higher-velocity shells 
(more darker points). It also shows that more packets are scattered from 
an initially red wavelength to blue and UV wavelengths (dark points with 
outgoing wavelengths < 3500 A below the diagonal). 
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Figure 5. Comparison of the integrated flux in different bandpass filters for the model series of SN 2001eh (left) and SN 2001ep (right). The curves show the 
change in magnitudes in a specific filter band (cf. Fig. [4) relative to the base model versus the change in composition in the outer shell of the model. 



to be scattered out of their original wavelength. In the base model 
10% of all packets emerging with A < 3750 A encounter their last 
interaction in a shell with v > 20 000 km s~ ; in the model with high 
Fe abundance this fraction is 20%. Counting all packets emerging 
from the envelope regardless of wavelength, the fractions of pack- 
ets encountering the last interaction at v > 20 000 km s~' are 3.9% 
(base model) and 5.6% (Fe enhanced). 

The condition for the reverse fluorescence process to lead to 
an increase of the UV flux is that the UV lines are saturated such 
that photons can escape from the lines without being re-absorbed 
immediately. Part of the UV enhancement seen in the models with 
increased Fe abundance is therefore caused by the presence of a 
larger number of saturated Fe lines in the UV. The decrease of UV 
flux in the models with enhanced 56 Ni abundance indicates that the 
dominating Co lines in those models are on average less saturated. 
Therefore the optical depth in those lines increases with abundance. 

The second consequence of the abundance variation in the 
outer shells is a change in the temperature and ionization structure 



caused by the change in line opacity. Fe group elements have a large 
number of lines in the UV and blue part of the spectrum. The high 
velocities cause them to overlap and effectively prevent photons to 
emerge. The higher radiation energy density in the ejecta leads to 
higher temperatures and consequently a higher state of ionization. 

Fig. [7] shows the mean radiation temperature (upper panel) 
and the ratio of doubly to singly ionized Fe (lower panel) for 
the SN 2001ep-based model series in which 56 Ni is varied above 
14 500 km s~' (shaded area). A change to lower abundances does 
not have a strong effect but with increasing 56 Ni abundance a 
clear shift to higher temperatures and consequently doubly ionized 
species is induced. Backwarming extends the region affected to ve- 
locities lower than j; s heii above which the actual variation has been 
imposed. 

IPinto & Eastman! d200Gt) discuss that in a gas having the phys- 
ical conditions typically found in supernovae (low matter densi- 
ties, high radiation energy densities) and dominated by sufficiently 
complex ions with a large number of possible line transitions, an 
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Figure 7. The upper panel shows the temperature as a function of velocity 
in the SN 200 lep model series where the abundance of 56 Ni was varied in 
the outer region above « s hcii = 14 500 km s (shaded area). The lower panel 
shows the relative ionization of Fe III/Fe II for the same models. 

The models differ slightly even below because of the different number 
of photons scattered back into the envelope from outer layers depending on 
the heavy-element abundance. 
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Figure 8. Sum of rs b for all lines in a wavelength interval of 10 A for the 
ions indicated in the plot legend. The Sobolev optical depths are derived for 
the model of SN 2001ep. 



equilibrium of the radiation field is established solely by radiative 
transitions. They regard this as an additional thermalization pro- 
cess that is not related to collisional and true continuum processes. 
In contrast to real thermalization, however, this equilibrium state 
does not couple the thermal properties of the gas to the radiation 
field. Therefore, the emitted spectrum depends only weakly on the 
temperature of the gas (as long as the ionization does not change 
significantly with temperature). Instead, the spectrum reflects the 
distribution of line strengths as a function of wavelength. A change 
in the dominant species in the outer layers, either by a variation of 
the composition or the ionization state of the gas, can affect the den- 
sity of spectral lines in the UV and may have a strong effect on the 
emergent spectrum. Fig.[8]shows the distribution of line strengths in 
terms of Sobolev optical depths for some of the relevant ions found 
in the mixture of SN ifl The qualitative difference of the distribu- 
tion between different elements explains why different wavelength 
regions are affected when Ti and Cr are varied compared to the 
variation of Fe and Co. Overall the higher ionization stages tend to 
have fewer strong lines in the relevant spectral region. Therefore, 
a higher ionization decreases the total opacity. This introduces a 
feedback effect that prevents the opacity from growing arbitrarily 
with the abundance of Fe-group elements. 

4.2 Variation of the density structure 

4.2.1 Setup 

For this series we kept the inner structure of W7 out to a velocity 
v shell = 15 000 km s~' and replaced the outer part with a power- 
law density p oc (u/u s h e iiy with power-law index /3 ranging from 
-5 to -16. The grey shaded area in Fig. [9] shows the range of 
variation considered. For compar ison, the density stru ctures of the 
pure ID deflagration model W7 dNomoto et al.lll984|) and the de- 
layed detona tion explosion m odels DD4 dWooslev & WeaveJl994h 
and WDD2 dlwamoto etalll 19991) are also plotted. The DDT mod- 
els have somewhat more mass at high velocities (i.e. a shallower 
density gradient). The variation considered in the parametrized 
models covers the typical range of gradients for deflagration and 
delayed detonation models. All models are extrapolated to c ~ 
70 000 km s - ' assuming power-law densities to ensure that the outer 
boundary does not affect the results. 

Results the f rom multi-dimensional explosion models of 
iRopke et alj J2007I) indicate that the "bump" and the relatively 
sharp change in the gradient of the W7 deflagration model around 
13 700km s _1 may be an artefact of the one-dimensional treatment. 
Multi-dimensional models tend to have a smoother density struc- 
ture, although this can depend strongly on the assumptions made 
when averaging multi-dimensional models over angles to obtain a 
ID density structure. If the deflagration phase is followed by a de- 
layed detonation, material is burned to higher velocities than in a 
pure deflagration because the flame is accelerated and can reach the 
outer regions of the white dwarf before the de nsities are too low for 
burning to proceed ( Ropke & Niemeyer 20q3). 

At the epoch considered here, a modification of the density 
gradient at high velocities is expected to affect the UV part of the 
spectrum most strongly because the region where the optical spec- 
trum is formed lies at much lower velocities. We restrict this study 
to one-dimensional explosion models because our radiative transfer 
models are spherically symmetric. 

7 The data shown in this plot are extracted from the SN 2001ep model but 
the distributions depend only weakly on the model. 
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Figure 9. The model setup for the variation of the density structure in the 
outer part. The W7 density (solid line) is replaced by a power-law with dif- 
ferent power-law indices p between —5 and —16 (shade d area). For com par- 
ison th e two delayed-detonation models WDD2 l lwamoto et al. 199jJ) and 
DD4 jWooslev & W eaver 1994) are also shown. For the spectral models the 
density structures are extrapolated to higher velocities (~70 000kms~') to 
avoid artefacts resulting from the outer boundary. 
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4.2.2 Results 

Fig. [10] shows the model series for different gradients as well as 
the model based on the W7 density structure. The region which is 
mostly affected by the change is the UV bluewards of about 4200 A. 
The rest of the spectrum remains largely unchanged. 

The W7 model corresponds roughly to the model with /3 = 
-12. Models with a comparable or steeper gradient than the W7 
structure show similar features that tend to become narrower for 
steeper gradients. In the models with shallower gradients the fea- 
tures tend to become broader and shift bluewards for a power- 
law index of p < -9. In even shallower models the structures 
become progressively less distinct and the flux in the UV below 
2750A drops significantly. The overall effect is similar in both 
model series except for the double peaked emission feature cen- 
tred around 3000 A, which shows a stronger variation in the series 
of SN 2001eh. 

Fig. II llshows the integrated flux in the bandpass filters shown 
in Fig.|4]for the model series with density variation (Fig.llOt. The 
strongest variation is observed for the three bluest filters (fl70w, 
fl85w, and f218w), while other bands show only minor changes. 
For very steep gradients the flux is enhanced in those bands rela- 
tive to the base model. As the gradient flattens the flux drops by 
up to 2.4mag in the fl85w-band. As expected, the results for the 
two supernovae are very similar. The difference at ~ 3000 A is not 
visible in the integrated flux because the f300w filter is too broad. 

We also computed models for the density structures of 
the two one-dimens ional dela yed detonation models W DD2 
dlwamoto et afl ll998) and DD4 dWooslev & Weavej|l994h using 
the same composition as in the base model for SN 2001eh. The 
resulting spectra are shown in Fig. [12] Both of the delayed detona- 



Figure 10. The full UV and optical spectra of the model series with varying 
power-law density above f s hell = 15 000 km s . The upper panel shows the 
models based on the model for SN 2001 eh, the lower panel shows those 
based on SN 2001ep. The colour coding refers to the power-law index fl 
used for each model. For comparison, the black line refers to the original 
W7-based model spectrum, which is close to the /? = —12 case in the outer 
part. The insert shows the UV part of the spectrum in logarithmic scaling. 

tion explosion models have a somewhat higher kinetic energy than 
W7 and thus more mass at higher velocities (cf. Fig. [9}. The den- 
sity gradients in the outer parts of both models are comparable and, 
therefore, the model spectra are not very different. Compared to the 
W7 model, the prominent features in the blue and UV are broader 
and shifted bluewards in the delayed detonation models. The same 
comparison for SN 200 lep shows qualitatively similar features but 
the difference to W7 is smaller. 

4.2.3 Discussion 

Line opacities scale with density. Therefore the modification of 
the density gradient in the outer layers affects the wavelength 
range over which a photon can interact with a spectral line. Shal- 
lower density structures increase the effective range of Doppler- 
shifts from the line-centre letting more lines overlap and leading 
to broader line features. This makes the line blocking at UV wave- 
lengths even more efficient. Fig. Q~3] shows the temperatures and 
the ionization ratio Fe III/Fe II for the SN 2001ep-model series 
with a varying density gradient. While the temperature varies only 
mildly for models with a density gradient steeper than W7, it be- 
comes much higher in the outer parts for the models with very shal- 
low density gradients. As in the models with varying composition 
this behaviour reflects the increased heating of the ejecta by the 
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Figure 11. Variation of the integrated flux in the filters shown in Fig. 1101 
relative to the original W7-based models. 
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Figure 12. Spectra of the one-dimensional delayed detonation explosion 
models WDD2 and DD4 as well as the deflagration model W7 assuming 
an identical composition for all models. The series shown is based on the 
model for SN 2001 eh. 



higher radiation density caused by the more efficient line blocking. 
The ionization pattern shown in the lower panel of Fig. [131 roughly 
follows the temperature, with the exception that models with very 
steep gradients show a significant shift to higher ionization stages 
in the outer parts beyond v ~ 17 000 km s -1 . Here the low densities 
become the governing factor that drives the ionization towards the 
doubly ionized species although the densities are too low to leave 
a significant imprint in the absorption features in the optical spec- 
trum. The pattern is very similar for the other elements. The role of 
the distribution of line strength discussed in the previous subsection 
applies to this model series in the same way. The varying spectral 
appearance in the UV can therefore be qualitatively understood by 
the shift in ionization for the elements that provide the major source 
of opacity. 
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v [km/s] 



Figure 13. Similar plot as in Fig.[7]for the models with varying density gra- 
dients based on the model of SN 200 lep. The colour coding of the models 
is the same as in Fig. 1 10| The W7-model roughly corresponds to fl = -12, 
however is somewhat shallower at lower velocities. 



5 CONCLUSIONS 

We analysed the spectra of two SNe la that have spectral cover- 
age in the UV wavelength bands. Although the spectra were taken 
at very similar epochs, their appearances are different. The models 
show that the typical temperatures in the ejecta of SN 2001eh are 
somewhat higher than in SN 200 lep. The shape of the spectrum of 
SN 2001 eh is determined by absorptions of doubly ionized species 
(Fe III and Co III) while the shape of the spectrum of SN 200 lep 
is determined by lines of singly ionized iron-group elements (Fe II 
and Co II). The difference is especially visible in the UV bluewards 
of Ca II H&K. Based on these spectral models we investigated the 
dependence of the UV flux on the metallicity and the density gra- 
dient in the outer layers of the ejecta. In one series of models we 
varied the abundance of groups of heavy elements; in another series 
we modified the density gradient in the outer layers and replaced 
the adopted explosion model with a power-law density structure 
that has different power-law indices. 

The model series with varying composition shows that under 
certain conditions the UV flux may actually increase with increas- 
ing abundance of heavy elements while other models show the op- 
posite behaviour. The physical process that dominates the forma- 
tion of the UV flux in SNe la is a reverse-fluorescence process in 
which photons are absorbed at long wavelengths and re-emitted at 
a shorter wavelength in the blue or UV. Therefore, the shape of 
the spectrum is determined by the distribution of spectral lines as a 
function of wavelength. This distribution varies slightly for differ- 
ent elements, but is very sensitive to different ionization stages. A 
temperature change that results in a shift of the ionization balance 
directly affects the efficiency of line blocking in the UV, which cou- 
ples back to the temperature. This feedback results in a non-linear 
response of the UV flux to changes in the metallicity. This wave- 
length region is therefore a valuable probe for the outer regions of 
explosion models, which do not leave a strong signature at longer 
wavelengths. To further study this effect quantitatively, a more re- 
alistic treatment of the temperature structure, consistent with the 
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non-LTE state of the gas, should be used. In this work we can only 
give a qualitative assessment of potential variations. 

The density gradient mostly affects the breadth of observed 
line features. Shallow density gradients increase the velocity region 
over which a line can efficiently absorb and re-emit photons. For 
very shallow density structures the densely clustered lines from Fe 
group elements in the UV cause individual absorption features to 
merge together resulting in a strong depletion of the flux. Steeper 
density gradients lead to narrower and more pronounced absorption 
and re-emission features. 

Overall, the model series of both supernovae show the same 
qualitative changes. The main differences seen in the UV affect 
the shape and strength of the re-emission peak at ~3500 A, blue 
of Ca II H&K. 

To conclude, a variation in composition and density of the 
high-velocity layers of the ejecta, which may be introduced by dif- 
ferent progenitor environments or different degrees of mixing in the 
actual explosion, can have a significant impact on the observed UV 
flux even at later epochs when the optical spectrum is formed at 
much deeper layers. Depending on the exact physical conditions, 
a variation of metallicity can either reduce or increase the UV flux 
without leaving a strong imprint on the spectrum at longer wave- 
lengths. For a better estimate of the magnitude and the scatter of 
the variation found in real supernovae, more UV observations of 
SNe la, preferably over a range of epochs, are needed. With the up- 
coming repair of STIS and the installation of the Cosmic Origins 
Spectrograph (COS) on HST, it should be possible to obtain addi- 
tional UV spectra of SNe la that would be critical for such studies. 
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